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LONG-TERM  GOALS 

The  long-term  goal  of  this  project  is  to  develop  and  apply  practical  methods  for  data  assimilation  to 
improve  the  short-range  prediction  of  mesoscale  ocean  variability. 

OBJECTIVES 

The  primary  objective  of  this  work  is  to  develop  an  ocean  data  assimilation  system  that  exploits  the 
strengths  of  both  the  ensemble-based  (e.g.,  Evensen  2003;  Houtekamer  and  Mitchell  1998;  Tippett  et 
al.  2003)  and  variational  (e.g.,  Bennett  2002)  approaches  to  data  assimilation.  The  first  step  in  this 
project  is  to  perform  a  comprehensive  inter-comparison  of  an  ensemble-based  data  assimilation  system 
with  a  4d-var  system  for  a  suite  of  coastal  model  configurations.  The  second  step  is  to  identify  the 
strengths  and  weaknesses  of  each  system  and  to  improve  both  systems  by  borrowing  components  from 
the  other  system.  Ultimately,  a  single  ensemble-var  system  will  be  developed.  We  will  investigate  the 
extent  to  which  the  ensemble-var  system  can  outperform  both  the  ensemble-based  and  variational 
approaches,  both  in  terms  of  forecast  skill  (accuracy)  and  computational  efficiency  (throughput). 

APPROACH 

An  ensemble-based  data  assimilation  system,  based  on  an  Ensemble  Optimal  Interpolation  (EnOI) 
scheme  (Oke  et  al.  2002;  Evensen  2003),  has  been  developed  and  tested  for  a  range  of  applications  by 
Dr.  Oke  (Oke  et  al.  2005;  2008;  2009;  2010;  Yin  et  al.  2010;  O’Kane  et  al.  2010),  and  a  four¬ 
dimensional  variational  (4d-var)  system,  based  on  the  representer  method  (Chua  and  Bennett  2001; 
Bennett  2002),  has  been  developed  and  tested  for  coastal  ocean  applications  by  Dr.  Kurapov  (Kurapov 
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et  al.  2007;  2009  2010). 


Conceptually,  both  the  ensemble  and  variational  approaches  perform  analogous  tasks.  Both  methods 
implicitly  generate  estimates  of  the  system’s  background  error  covariance,  both  interpolate  and 
extrapolate  background  innovations  (model-observation  differences)  onto  the  full  model  state 
(including  all  variables  at  all  model  grid  points),  and  both  seek  to  minimize  some  norm  of  the  model- 
observation  misfits.  The  ensemble-approach  has  an  advantage  that  it  uses  the  non-linear  model 
operator  to  generate  and  evolve  the  system’s  error  covariance,  but  has  the  disadvantage  that  it  is  rank 
deficient  (i.e.,  usually  uses  too  few  ensemble  members  to  completely  span  the  error  sub-space  of  the 
system).  The  variational  representer-based  approach,  by  contrast,  is  full  rank,  but  it  uses  linear 
operators  to  derive  the  system’s  error  covariance,  or  representers.  The  ensemble-approach  is  typically 
implemented  as  a  three-dimensional  interpolation  (e.g.,  Oke  et  al.  2005;  2008)  and  the  representer- 
approach  is  fundamentally  four-dimensional  (Kurapov  et  al.  2007,  2009).  There  are  pros  and  cons  of 
both  approaches.  Under  this  project,  we  will  exploit  the  strengths  of  both  approaches,  with  the  goal  of 
developing  a  superior  ensemble-var  data  assimilation  system. 

A  key  hypothesis  in  this  study  is  that  the  adjoint  and  tangent  linear  model  developed  for  a  particular 
ocean  general  circulation  model  -  in  this  case  the  Regional  Ocean  Modelling  System  (ROMS)  -  can 
successfully  be  applied  to  a  different  ocean  model  -  in  this  case  the  Sparse  Hydrodynamic  Ocean  Code 
(SHOC;  Herzfeld  2006).  If  this  can  be  demonstrated,  the  technical  developments  of  the  ocean 
modelling  community  can  potentially  be  consolidated  significantly  -  and  the  progress  of  this 
community  could  be  accelerated.  We  have  made  steps  towards  testing  this  hypothesis  by  configuring 
both  SHOC  and  ROMS  for  the  same  region,  on  the  same  grid. 

The  EnOI  system  to  be  used  under  this  project  is  routinely  applied  to  regional  configurations  of  SHOC. 
The  technical  aspects  of  applying  the  4dVar  system  to  SHOC  have  been  achieved.  This  has  involved 
porting  the  Advanced  Variational  Regional  Ocean  Representer  Analyzer  (AVRORA;  Kurapov  et  al., 
2007;  2009;  2020)  to  the  same  computer  platform  as  SHOC,  and  enabling  AVRORA  to  assimilate  data 
from  the  same  database  that  the  EnOI  system  draws  on.  These  steps  have  been  largely  technical  -  but 
have  enabled  the  Australian-based  researchers  involved  in  this  project  to  acquire  a  good  working 
knowledge  of  AVRORA.  We  are  now  well-placed  to  fulfill  the  objectives  of  this  project  by  first  inter¬ 
comparing  EnOI  with  4dVar,  then  bringing  together. 

A  key  characteristic  of  the  4dVar  approach  is  the  state-dependence  of  the  background  error  covariance 
estimates  -  quantified  by  the  representers.  Enabling  state-dependent  estimates  of  the  background  error 
covariance  using  the  EnOI  is  one  way  of  bringing  these  approaches  together. 

Dr.  Peter  Oke  is  the  P.I.  on  this  project  and  leads  the  ocean  data  assimilation  activities  at  CSIRO 
Marine  and  Atmospheric  Research.  Dr.  Oke  is  working  closely  with  Dr.  Kurapov,  from  OSU,  under 
this  project,  and  also  in  collaboration  with  Dr.  Terrence  O’ Kane  and  Dr.  Chaojiao  Sun  from  CSIRO. 

WORK  COMPLETED 

A  series  of  comparisons  between  ROMS  and  SHOC,  configured  for  the  Bonney  coast  off  southern 
Australia  have  been  conducted.  Results  show  that  both  models  produce  comparable  wind-driven 
circulation.  The  experiments  include  both  ROMS  and  SHOC  nested  within  fields  from  a  global  ocean 
reanalysis,  and  fields  from  a  regional  EnOI-based  analysis.  Comparisons  demonstrate  that  both  ROMS 
and  SHOC  produce  more  realistic  circulation  when  they  are  nested  within  the  EnOI-based  analyses. 


The  next  step  of  applying  the  4dVar  system  to  both  models  is  underway.  Modules  have  been  written  to 
allow  the  4dVar  system  to  assimilate  observations  from  the  database  used  by  the  EnOI  system.  This 
will  facilitate  a  direct  inter-comparison  of  methods. 

State-dependent  estimates  of  the  background  error  for  the  EnOI  system  have  been  produced  by 
establishing  a  simple  EnOI-breeding  capability  (O’ Kane  et  al.  2010).  The  EnOI-breeding  approach  has 
been  applied  to  a  regional  model  of  the  Tasman  Sea.  Results  demonstrate  that  this  efficient  approach 
yields  state-dependent  estimates  of  the  system's  background  errors  that  are  meaningful. 

RESULTS 

A  comparison  between  ROMS  and  SHOC,  both  configured  for  the  Bonney  coast  off  southern 
Australia,  has  been  performed.  The  results  show  that  both  models  reproduce  the  variability  associated 
with  wind-driven  upwelling  (Figure  1).  The  intensity  of  the  upwelling  is  similar  in  both  models.  Both 
models  benefit  from  being  nested  within  the  EnOI-based  regional  analysis  (labelled  BODAS  BCs).  In 
these  experiments  the  EnOI-based  regional  analyses  (BODAS  BCs)  are  computed  by  assimilating 
observed  sea-surface  temperature,  along-track  sea-level  anomalies,  and  in  situ  temperature  and  salinity 
from  Argo  profiles  every  day.  The  background  field  for  each  analysis  is  derived  from  version  2.2  of  the 
Bluelink  ReANAlysis  (BRAN2.2;  Oke  and  Griffin  2010). 
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Figure  1:  Standard  deviation  of  sea  surface  temperature  for  February  2008 from  SHOC  with  BRAN 
boundary  conditions  (left),  ROMS  with  BRAN  boundary  conditions  (middle),  and  ROMS  with 
BODAS(i.e.,  regional  EnOI-based  analyses)  boundary  conditions. 

[Figure  showing  the  similarity  between  two  different  ocean  models  off  the  Bonney  coast,  and  the 
benefit  of  assimilated  within  EnOI-based  analysis  fields.] 

Breeding  is  a  simple  method  that  involves  the  perturbations  of  an  analysis  field,  and  the  sequential 
rescaling  of  the  perturbed  forecast  anomalies  with  respect  to  a  specified  norm.  In  theory,  breeding 
generates  the  perturbations  fields  that  project  most  strongly  onto  the  fasting  growing  modes.  Bred 
vectors  can  therefore  provide  an  estimate  of  the  likely  growth  of  the  background  error  over  a  given 
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assimilation  cycle.  The  background  errors  are  state-dependent.  An  example  of  breeding,  applied 
together  with  an  EnOI  system  (Oke  et  al.  2008),  is  presented  in  Figure  2  (from  O’Kane  et  al.  2010). 
Figure  2  shows  the  mean  of  4-bred  vectors  for  three  different  times,  along  with  an  estimate  of  the 
background  error  for  sea-level  (computed  by  differencing  a  forecast  with  a  verifying  analysis).  In  each 
case  shown  in  Figure  2  the  bred  vectors  have  large  amplitude  in  regions  that  correspond  to  large 
background  errors  (the  mean  of  the  bred  vectors  is  contoured  in  black  over  the  estimates  of  forecast 
error).  Note  that  the  bred-vectors  are  computed  without  prior  knowledge  of  the  actual  background 
errors.  The  growth  of  each  bred  vector  is  therefore  due  to  growing  instabilities  in  the  model  that 
correspond  to  the  development  or  evolution  of  mesoscale  eddies  and  the  associated  sub-mesoscale 
circulation.  In  practice,  such  instabilities  are  poorly  predicted  owing  to  their  chaotic  nature,  so  the 
background  errors  around  these  growing  instabilities  are  often  relatively  large. 
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Figure  2:  Examples  of  (a-d)  forecast  error  for  sea-level  ( colour )  and  the  ensemble  averaged  bred 
vector  (contours,  showing  plus/minus  0.35  m  contours )  in  the  Tasman  Sea;  and  (e-h)  four  bred 
vectors  overlaid.  Each  bred  vector  is  a  different  colour. 

[Figure  showing  the  alignment  of  the  bred  vectors  and  the  state-dependence  of  background  errors 

from  a  data  assimilating  ocean  model.] 


IMPACT/APPLICATIONS 

We  are  a  step  closer  to  evaluating  the  hypothesis  that  an  adjoint  and  tangent-linear  model  from  one 
model  (in  this  case,  ROMS)  can  be  used,  together  with  a  4dVar  system,  to  assimilate  observations  into 
a  second  independent  model  (in  this  case,  SHOC).  If  this  hypothesis  can  be  verified,  the  technical 
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challenges  of  developing  an  adjoint  model  and  a  tangent  linear  model  for  every  new  version  of  an 
ocean  model  will  be  significantly  simplified. 

The  breeding  method  employed  here  (see  Figure  2)  reliably  predicts  the  state-dependence  of  the 
system's  background  error.  When  combined  with  EnOI,  or  a  4dVar  system,  breeding  should  ultimately 
make  these  systems  more  efficient.  Breeding  is  also  suitable  for  adaptive  sampling,  where  additional 
observations  (e.g.,  gliders)  could  be  deployed  in  regions  of  large  errors  or  growing  instabilities. 

RELATED  PROJECTS 

Bluelink  is  a  partnership  between  CSIRO,  the  Bureau  of  Meteorology  and  the  Royal  Australian  Navy. 
Many  of  the  research  activities  undertaken  in  Bluelink  have  strong  synergies  with  the  project  that  is  the 
subject  of  this  annual  report.  The  main  objective  of  Bluelink  is  the  development  and  application  of  an 
ocean  forecast  system  for  the  mesoscale  circulation  around  Australia.  Applications  of  the  Bluelink 
system  are  well  documented  (e.g.,  Oke  et  al.  2005;  2008;  2009;  2010;  Schiller  et  al.  2008).  An 
indication  of  the  performance  of  BRAN  is  shown  in  Figure  3,  where  observed  and  modelled  sea-level 
is  compared  for  several  unassimilated  tide  gauges  off  south-eastern  Australia  (Oke  and  Griffin  2010). 
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Figure  1:  Time  series  of  (A)  meridional  (bold)  and  zonal  (thin)  wind  stress  at  35. IS;  and  (B-F)  sea- 
level  anomaly  from  BRAN  (thin)  and  unassimilated  coastal  tide  gauges  off  south-eastern  Australia. 
Tide  gauge  locations  are  denoted  by  the  dots  in  the  panel  to  the  right. 

[Figure  showing  the  good  agreement  between  observed  and  modelled  sea-level  off  eastern  Australia 

from  the  Bluelink  Reanalysis  system.] 


The  most  recent  activities  under  Bluelink  involve  the  development  of  a  near-global  (excluding  the 

Arctic)  eddy-resolving  ocean  model.  The  new  Bluelink  model  will  underpin  future  ocean  reanalyses 

and  Australia's  future  operational  short-range  ocean  forecast  system.  Application  of  the  EnOI  system 

referred  to  throughout  this  annual  report  is  underway. 
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